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Abstract 

A variety of complex fluids under shear exhibit complex spatio-temporal behaviour, including 
what is now termed rheological chaos, at moderate values of the shear rate. Such chaos associated 
with rheological response occurs in regimes where the Reynolds number is very small. It must thus 
arise as a consequence of the coupling of the flow to internal structural variables describing the 
local state of the fluid. We propose a coupled map lattice (CML) model for such complex spatio- 
temporal behaviour in a passively sheared nematic liquid crystal, using local maps constructed so 
as to accurately describe the spatially homogeneous case. Such local maps are coupled diffusively 
to nearest and next nearest neighbours to mimic the effects of spatial gradients in the underlying 
equations of motion. We investigate the dynamical steady states obtained as parameters in the 
map and the strength of the spatial coupling are varied, studying local temporal properties at a 
single site as well as spatio-temporal features of the extended system. Our methods reproduce the 
full range of spatio-temporal behaviour seen in earlier one-dimensional studies based on partial 
differential equations. We report results for both the one and two-dimensional cases, showing that 
spatial coupling favours uniform or periodically time-varying states, as intuitively expected. We 
demonstrate and characterize regimes of spatio-temporal intermittency out of which chaos develops. 
Our work suggests that such simplified lattice representations of the spatio-temporal dynamics of 
complex fluids under shear may provide useful insights as well as fast and numerically tractable 
alternatives to continuum representations. 

PACS numbers: 05.45.Ra,05.45.Pq,83.30.Xz 



* Electronic address: 'kaniil@ims c.res . iiil 
^Electronic address: menon@imsc.res.inl 
■f Electronic address: lsudeshna@imsc . res .inl 



I. INTRODUCTION 



Unusual dynamical steady states are obtained in a large number of experiments on com- 
plex fluids driven out of equilibrium [l|-|6|. When such fluids are sheared uniformly, the shear 
stress 0" is typically regular at very small shear rates 7. However, at larger shear rates the 
response is often unsteady, exhibiting oscillations in space and time as a prelude to inter- 
mittency and chaos |7Hll|. In this non-linear regime, complex fluids under shear exhibit a 
variety of instabilities, including instabihties to "shear banded" states|12Nl7|. Such banded 
states arise from an underlying multi-valued constitutive relation connecting the stress and 
the shear rate, and are often obtained as a precursor to spatio-temporal intermittency and 
chaotic behaviour in flow response 18|-|26|. 

Such rheological chaos must be a consequence of constitutive non-linearities, since 
Reynolds numbers associated with the flow are too small for convective non-linearities to 
be important 16|, ll7|]. Such constitutive non-linearities originate in the non-trivial internal 
structure of the fluid and its coupling to the flow. Recent rheological studies of "living poly- 



mers o 
value [7 



3tain an oscillatory stress response to steady shear at shear rates above a 



10|. Such an oscillatory response turns chaotic at still larger shear rates 



threshold 

7Hio|. It 



has been argued that a hydrodynamic description of this behaviour requires coupling the 
internal orientational state of such a polymeric fluid to the flow, motivating the study of the 
problem addressed in this paper 27|, |28|. This is the model problem of the spatio-temporal 



description o 



shear flow 
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an orientable fluid, such as a nematic liquid crystal, placed in a simple steady 



3l|. 



There is a substantial body of previous work on the dynamical states of complex fluids 
under shear. A model due to Fielding and Olmsted expresses the stress as a function of 
a microstructural parameter chosen, for illustrative purposes, to be the micellar length, 
which itself evolves in response to the shear rate. The microstructural parameter yields 
a viscoelastic contribution to the stress, over and above the regular fluid contribution 32 1 . 
Fielding and Olmsted show that their model exhibits spatio-temporal rheochaos. Aradian 
and Gates have proposed a one-dimensional model for the instabilities of a shear-banding 
fluid system, writing down an equation for the time-variation of the shear stress which 
depends both on the instantaneous value of the strain rate as well as on the previous history 
of the stress|3S]. This single non-local equation can be cast as two coupled local equations. 



one for the stress as well as another for a "memory" term, arising out of the single equation 
for the stress evolution. This simple model yields regimes of periodic as well as chaotic 
behaviour J3 31] ■ 



Both these models assume simplified scalar descriptions of the internal microstructure. A 
recent, comprehensive study of a shear-banding interface by Fielding and Olmsted, based on 
the diffusive Johnson-Segalman (DJS) model, shows that the interaction of multiple shear 
bands can yield a time-dependent stress response possessing attributes of low- dimensional 
chaos 3J]. However, such approaches do not examine how such a stress response might 
arise from an underlying microscopic equation of motion. Recent work by Chakraborty, 
Dasgupta and Sood on a one-dimensional model for nematic rheochaos extends the model 



of Refs. 271, |28| by incorporating hydrodynamics, finding stable shear banding as well as the 



coexistence of banded and spatio-temporally chaotic states 35|. Further, the DJS model is 
derivable as a specific limit of their model, in which the equation for the order-parameter 
part of the stress is linearized about the isotropic limit. 

In this paper, we present results from a comprehensive study of a simple coupled map 
lattice model for rheological chaos, as appropriate to nematic systems under steady shear. 
Our local "microstructural" variable represents the orientation and degree of coarse-grained 



m 
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35|. We compute the 



order of nematic molecules in the flow, as in the work of Refs. 
contribution to shear stresses arising from the evolution of this local variable, showing how 
uniform, periodic and spatio-temporally chaotic behaviour in this quantity can be accessed. 
The use of coupled map lattices to represent, at a coarse-grained level, behavior of in- 
trinsically non-linear dynamical systems coupled in space is at least two decades old 37| • 
Coupled map lattices provide relatively simple models whenever it can be assumed that the 
dynamics can be naturally decoupled into a dominant local dynamics representing behaviour 
at a single point in space (or small coarse-grained region) and a spatial coupling term which 
connects this local dynamics weakly across spatial locations. The coupling term idealizes 
gradient terms in the underlying continuum equation of motion. Coupled map lattices are 
well-suited for computer simulations, since they are naturally discrete in space and time. 
(Experimental data are, in fact, close to the CML situation, since any real-life measurement 
requires discrete sampling of the underlying time evolution and every experiment has some 
minimum threshold for spatial discrimination, providing a lattice scale.) Coupled map lat- 
tices have been used with success by several authors in the study of phase-ordering problems 



as well as in a host of other applications |36l. |37] . 

We begin by constructing a local map for nematics under shear, obtained by discretizing 
a set of coupled ordinary differential equations (ODE's) describing the continuous time, 
spatially local version of this dynamics. These local equations have been shown to exhibit 
periodic and regular regimes as well as chaotic regimes. We benchmark this map through 
a detailed comparison to the results from the study of the ODE system, showing that 
the qualitative and quantitative aspects of the phase diagram in this single site limit are 
rendered accurately. We then generalize this to the spatially coupled case by connecting 
nearest neighbour maps in a specified manner. The shear enters at the level of the local 
map, where it is specified in terms of a single parameter. We take the point of view that the 
complexity of the spatio-temporal behavior in the physical problem can be captured by the 
most elementary version of spatial coupling, which, for simplicity and following virtually all 
work on coupled map lattices, we take to be diffusive |37l|. 

This local map is shown, in agreement with previous work, to exhibit a large number of 
complex phases, including uniform (flow aligning in the nematic), tumbling, kayaking and 
chaotic phases, in addition to phases which combine one or the other of these attributes |38l. 



391 ]. While the nematic responds to the fluid through flow alignment as well as reactive 
and dissipative terms in the equation of motion, we make the approximation of ignoring 
the back-reaction of changes in nematic order on the fluid. Thus, our approach omits the 
hydrodynamic interaction, since we assume that the flow always remains passive. This is 
a major assumption. However, it does have the virtue that a variety of spatio-temporal 
phenomena with relevance to both the experiments as well as to earlier modeling exercises 
can be demonstrated to exist in this simple system and are amenable to analysis. 

Our second approximation is that we study, for the most part, simple diffusive couplings 
between sites, ignoring the advective terms. Consistent with this, we use simple periodic 
boundary conditions on the local field. (We would otherwise have had to implement a more 
complex Lees-Edwards boundary condition on the fields and ensure an appropriate anchoring 
condition at the boundaries [40|].) Thus, in our model, the shear enters the local dynamics 
but its effects are ignored at larger scales. We pursue this line of investigation because 
our interest is specifically in the effects of including spatial couplings into a model which 
provides an accurate description of the temporal behavior of sheared nematics assuming 
spatial behavior to be uniform. We believe, and in some cases have tested this assumption. 



that incorporating the simplest form of spatial coupling should be sufficient for us to be able 
to explore the full spatio-temporal complexity of the sheared nematic problem. 

The outline of this paper is the following: Section II outlines our numerical methods for 
the construction of the local map. We begin by providing the local equation of motion for a 



passively sheared fluid of nematogens, following the work of Refs. [38|, |39| . To enforce sym- 
metry and tracelessness, it is customary to project these (tensor) equations onto a suitable 
tensor basis. We then construct, through a simple Euler discretization, a map within this 
basis, showing that it can be used to obtain all the states obtained by ODE-based methods 
for this problem. The following section. Section III, describes the construction of the cou- 
pled map lattice, illustrating how the local maps constructed in Section II can be coupled in 
space, in both one and two dimensions, along standard lines. Section IV describes our results 
in the one-dimensional case, examining the effects of spatial coupling in both regular and 
complex regions of the local map. Section V describes our results for the two-dimensional 
case, studying, as in the one-dimensional case, the behaviour in both regular and complex 
regimes of the local phase diagram. Section VI contains a discussion of our results as they 
relate to a quantification of spatio-temporal complexity in our model, while Section VII 
contains the conclusions of this study. 

II. A LOCAL MAP FOR NEMATODYNAMICS 

We begin with the continuum equations of motion for a nematic in a specified fiow 
field. These equations use the tensor representation of the order parameter in a nematic. 
In thermal equilibrium, such order parameter configurations are weighted by a Landau- 
Ginzburg-de Gennes free energy. In a specific Cartesian tensor basis, these equations, in the 
approximation that spatial fiuctuations in nematic order are absent, can be cast in terms 
of equations of motion for five expansion coefficients, corresponding to the five independent 
parameters characterizing a real symmetric traceless tensor. These equations of motion, 



which are ordinary differential equations (ODE's), are recast as a map, as shown below 41 1. 



A. Equation of Motion for Nematics 

The derivation of the nonlinear relaxation equations for the symmetric, traceless second 
rank tensor Q characterizing local order in a sheared nematic is available in earlier work 
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31 



42 



49|. The order parameter is often conveniently expressed as 

Qa/3 = — ( naUp - -5af3 j + — {IJfB - ^^711^) , (l) 

where the director n is defined as the normalized eigenvector corresponding to the largest 
eigenvalue of Q, the subdirector 1 is associated with the sub-leading eigenvalue, and their 
mutual normal m is obtained from n x I. The quantities Si and S2 represent the strength 
of uniaxial and biaxial ordering: |si| ^ , S2 = is the uniaxial nematic whereas Si, S2 7^ 
with S2 < 3si defines the biaxial case 5G|. 

Defining b := ^(b + b"^) — ^(trh)6 to be the symmetric-trace 



ess part of the second-rank 



tensor b, the equation of motion for Q in a velocity field is [29|, 

^ - 2frQ - 2a'f^ + ro^$ = -v^^T, (2) 

at '^ Ta 

where the tensor ft = ^((Vv)-*^ — Vv), T = |((Vv)-^ + Vv) and Vv is the velocity gradient 
tensor, with v = 'jye^, where e'' is a unit vector in the x— direction. The velocity is 
along the x direction, the velocity gradient is along the y direction, while z is the vorticity 
direction. The quantities r^ > and Tap are phenomenological relaxation times, a' describes 
the change of alignment caused by T and $ = d(j)/dQ,, with the free energy 0(Q) given by 

0(Q) = IaQ : Q - ^-VGBiQ ■ Q) : Q + ^C(Q : Q)^. (3) 

If the spatial variation is also taken into account, V^Q and VV ■ Q should also be included in 
the above expression . The notation Q : Q represents QijQji, with repeated indices summed 
over. Here A = Aq{1 — T* /T), and B and C are constrained by the conditions Aq > 0, 
fi > 0, C > and 52 > |AoC. 

The symmetric traceless alignment tensor Q has five independent components. Assuming 
spatial uniformity, so that gradients of the Q tensor can be dropped, a system of 5 coupled 
ordinary differential equations (ODEs) for the 5 independent components of Q can be ob- 
tained with the choice of a suitable tensor basis. Choosing the standard orthonormalized 
Cartesian tensor basis leads to the expansion 



Q = ^ a,T^ (4) 



A;=0 



with 

T^ = v^e^J/, T^ = V2e^^ T^ = V2e^^. (5) 

B. Dynamics of Sheared Nematics from a Local Map 

We work in the tensor basis described above, representing the equations of motion of Eq. [2] 
in terms of the coupled equations of motion for the five coefficients Oq . . . 04. The problem of 
representing the time updates in terms of a local map is most easily approached by consid- 
ering the lowest order Euler discretization of the underlying differential equations. (There 
are alternative methods of constructing maps from local dynamics governed by ODE's, in- 
cluding stroboscopic methods and methods which use Poincare sections; however, the choice 
we have made is the simplest given the variety and complexity of the dynamical states we 
would like to describe.) 



Scaling parameters as in Ref. 38|, [391], and making the same choice of numerical values 



as in Ref. 



27 



28| . we obtain the following map 



3 



/o(a*) = ao*+ A ( -(2a' - 3ao)ao - 3(af + a^) + -(03 + K 
/i(ai*) = ai*+ A ( -(2a^ + 6ao)ai + 702 + - V3(a^ - a^) j 

/2(a2*) = (32*+ A ( -(2a^ + 6ao)a2 - 701 + 3A/3a3a4 + —\ki 



a\az + 0213-4) 



hW) = a3*+ A ( - {2a^ - 3ao)as + ^7«4 + SV3{ 

/4(a.4*) = 0-1+ A - (2a' - 3ao)a4 - -703 + 3A/3(a2a3 - 0104) . (6) 

Here t indicates discrete time steps, and {.}* denotes the value of the quantity {.} at time 
step t. All the functions / denote the locally updated value of their argument at a time step 
(t -|- 1). For the purely local map, /i(a*) = o*^^; however, for the coupled map, the value of 
/i(a*) is computed as an intermediate step, prior to the diffusive step which yields the final 
quantity a*^"^. 

The quantity a = Oq + a' + a' + a' + a'. We choose A = 0.01 for all our calculations, but 
have checked that changing A by upto an order of magnitude does not affect our results. 

8 



Our choice of parameters implies that the system in the absence of shear is at the hmit of 
metastabihty of the isotropic phase. Our choice of the value for A captures all the features 
of the full local phase diagram obtained in Refs. J38|, |39 |. 

The order parameter part of the stress is proportional to contributions from the Landau- 
de Gennes free energy as well as from the gradient terms, which we represent through the 
spatial coupling term in the coupled map lattice. This is obtained as described in the 
following sections. 

C. Phase Behaviour of the Local Map 

Examining the dynamical steady states of this map at a large number of points in the 
space spanned by (7, A^) yields a complex phase diagram admitting many states - aligned, 
tumbling, wagging, kayak-wagging, kayak-tumbling and chaotic - as functions of the shear 
rate 7 and a phenomenological relaxation time A^ which is a parameter in the equations 



of motion 



38 



39 



51| . Fig. dl]) exhibits the dynamical states found in the map for the 
uncoupled case, in terms of a phase diagram in the quantities A^ and 7. Such a phase 
diagram bears considerable similarities to phase diagrams obtained by other authors in the 



PDE representation; see, for example. Fig. 7 of Ref. [i 

The states in this phase diagram are labeled as follows: The first is the state labeled 
A, which is the Aligned state, where all dynamics ceases, and the director is aligned at 
an angle to the flow. In the standard Couette geometry, the velocity field and the velocity 
gradient form a plane, called the vorticity plane. In our case, this is the x — y plane. If 
the director lies in the vorticity plane and rotates about an axis (the z— axis) perpendicular 
to this plane, the dynamical state is called a Tumbling state. The tumbling state is 
denoted by T in the phase diagram of Fig. ([T]). If the director, while lying in the plane, 
executes oscillations, the dynamical state is called a Wagging state. The wagging state is 
represented in the local phase diagram by the symbol W. 

If the director rotates and oscillates, moving out of the vorticity plane, the dynamical 
states are called Kayak- Tumbling and Kayak- Wagging respectively. They are repre- 
sented as KT and KW in the local phase diagram. If the dynamics is a mixture of complex 
intermittent behaviour and coexisting attractors, the state is called Complex and is repre- 
sented by C in the phase diagram. Clearly the interesting region in the phase diagram lies 
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Y 3.0^ 



FIG. 1: [Color Online] Phase diagram for dynamical behaviour in the local map defined through 
EqnsEl with the parameter Afc plotted on the x axis and 7 on the y axis. Here T denotes the 
tumbling state, W the wagging state, KT the kayak-tumbling, KW the kayak-wagging state, A 
the aligned state and C the state in which complex dynamics is seen. These states are discussed 
further in the text. 

in and near the region labeled C 

Fig. [1] is obtained in the following way. The phase-space of the 7 and A^ variables is 
gridded and an initial random initial condition chosen at each point. After the passage 
of an initial transient state, the system goes to dynamical attractors, ranging from simple 
spatiotemporal fixed points to complex intermittent behaviour. These dynamical attractors 
are identified with one of the states described above, i.e. A, K, T, KW, KT or C. In some 
regimes, one sees a coexistence of states i.e. KT and T and KW and W i.e. different initial 
conditions can give rise to different asymptotic behaviour in the long time limit. 

III. COUPLED MAPS FOR NEMATODYNAMICS 

Our spatially coupled model is built up from the local maps given in Eq. |6l These maps 
are placed on the sites of a regular lattice in one and two dimensions and can be coupled 
via several different coupling schemes, as described below. The generalization to arbitrary 
dimensions as well as different coupling schemes is a straightforward one. 

For a one dimensional lattice, with sites indexed by the label i, the five variables 
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(ao(i) . . . a4(i)) on each lattice site evolve in discrete time t as: 

</'^ = </': + ^e(0*+i + 0*_i-20*), (7) 

where G (a*+\ a*i+\ 4+\ 4+\ 4+^) and 0' G (/oK)> /i(«i), /2(a2), /sK), /4K)) -^eree 
is a coupling constant which is chosen to take values between and 3/2. 

For the two dimensional case we consider a square lattice with site index {i,j) and with 
the set of five variables {ao{i,j), ai{i,j), a2{i,j), a3{i,j), a4(z,j)) on each lattice point at 
time-step t evolving in time as : 

+— e (0-+1J+1 + 0--1J-1 + 0*-ij+i + 0U1J-1) - e0^' (8) 

where G (4+\ a*i+\ 4+\ a3+\ 4+^) and 0' G (/o(a*o), /i(a*i), /2(4), /3(4), ^(aD) , and e 
is a coupling constant having value between and 1. The choice of the numerical coefficients 
1/6 and 1/12 in the coefficients of the nearest and next-nearest neighbour terms are standard 
choices in the CML literature. They represent choices of lattice discretization which are as 
close as possible to the continuum limit. 

The local value of the shear stress {(Jxy)ij at the (two-dimensional) site {i,j) is obtained 



from the following definition 52|] : 



or. 



V2 
^xy)iJ +~^^ {(Pi+lJ + (t>i-lj + (Pi,j+l + 4>i,j-l) 



xyji,] \'xy/^,j ■ r. 



a/2 

+— e (</>i+i,j+i + 0i-i,j-i + 0i-i,i+i + 0i+ij-i) - V2e0M, (9) 

where {(Jxy)i,j is given by 

{a',y)ij = {2V2a2a^ - 6^6(0304/2 - (a2ao/v^)))i,i, (10) 

and (f) = a2. 

The definition for the one-dimensional case follows from: 

V2 

{(Txy)i = «j,)i +^e (<^i+l + (t>i~l - 20i)' (11) 

where {cr'^y)i is given by 

(a'^y), = {2V2a2a^ - 6^6(0304/2 - (0200/^8)))^, (12) 

11 



and = 02- 

We have also experimented with other choices of the update rule. While the update rule 
of Eq. [7] can be termed as the post-update rule, in which the terms on the right hand side 
are calculated at time t + 1, one could alternatively use the pre-update rule, in which these 
terms are calculated at time t. We have checked that varying this choice of update rule 
does not affect our results. In the equation for the two-dimensional update, (Eq. |8]), we 
have checked that dropping the next-nearest neighbour term also does not affect our results 
significantly. Thus, a variety of possible update schemes appear to yield consistent results 
for the spatio-temporal behaviour of our coupled map lattice, underlining the generic nature 
of our results. 

In the sections below we describe the relevant features emerging from the simulation of 
the coupled map lattice above for a wide range of parameters and initial states. 

IV. THE ONE-DIMENSIONAL COUPLED MAP LATTICE 

In this section, we describe our results for the one- dimensional case, concentrating on the 
effects of the inter-site coupling, both within and outside the regime labelled C in the phase 
diagram. 

A. Local dynamics 

Fig. (J2]) shows the dynamical phases exhibited by a generic site randomly chosen from 
the one dimensional ring. The sites are coupled according to the scheme given in Eq. |S1 with 
coupling constant e = 0.1 (a) and e = 0.5 (b). 

It is evident from comparisons with Fig. ([1]) that the local dynamics of a generic site in 
the coupled system is similar to the uncoupled case. This indicates that spatial coupling 
does not alter the nature of the local dynamics qualitatively. The most significant influence 
of spatial coupling occurs near the C region, which appear to be somewhat broadened 
with spatial coupling, while the coexistence regimes are reduced in size. In addition, the 
fairly uniform KT state is now "studded" with points displaying complex behaviour. This 
indicates coexistence of complex and KT behaviour, with certain initial states leading to 
complex dynamics, while others lead to a uniform KT state. (It is difficult to determine 
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FIG. 2: [Color Online] Phase diagram summarizing the local dynamical behaviour of the spatially 
coupled one-dimensional system, with A^ plotted on the x axis and 7 on the y axis. As before, T 
denotes the tumbling state, W the wagging state, KT the kayak-tumbling, KW the kayak-wagging 
state, A the aligned state and C the state of complex dynamics. The spatial coupling constant 
e = 0.1 (a) and 0.5 (b), for a ring of 200 lattice points. The location of the states in the phase 
diagram is largely similar to that of Fig. [T] with the exception that, at isolated points, mainly 
within the KT phase, one sees complex behaviour. The regime in the phase diagram occupied by 
the C phase shrinks as the coupling constant e is increased. 

whether the complex behaviour we see is a very long transient or true asymptotic behaviour.) 
The tumbling T and wagging W regions, however, are very stable. 

1 . Local behaviour of regular regions 

Figs. ([2])-(Ill) show the value of the scalar order parameter si, the biaxiality parameter 
S2 and the z component of the director n. Fig. ([3]) is obtained using parameter values 
appropriate to the T and W regions of the local phase diagram, with a coupling constant e = 
0.1. This displays completely regular behaviour, with these quantities varying periodically 
keeping the director in the vorticity plane. Fig. (jl]) is obtained using parameter values 
appropriate to the KT and KW regions of the local phase diagram and indicate that the 
director can now fluctuate out of plane whereas all quantities vary smoothly and periodically. 

Figure ([5]) shows the local time period with which these quantities fluctuate as a function 
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of 7 at Afe = 0.9. The time period is inversely proportional to shear rate, as is evident from 
the fit to the data points. 
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FIG. 3: [Color Online] Local dynamics in a ring of 200 lattice points, with coupling e = 0.1, showing 
the temporal evolution of si, ,32 and n^- These are displayed for (a) si ,S2 and n^ with (a) Xk = 0.9 
and 7 = 5.0 (b) Xk = 1-1 and 7 = 5.0, and (c) Xk = 0.9 and 7 = 4.0. These states are all drawn 
from the T and W parts of the local phase diagram. Note that 77,2 = in all these states whereas 
si and S2 are periodic functions of time. 



2. Local Dynamics in the Complex Region 

The local behaviour in the complex region, denoted by C in the local phase diagram, is 
exhibited in Fig. ([6]), which shows si and S2 and Uz- The results suggest that the sites display 
intermittent behaviour. These results are obtained for parameter values at the boundary of 
the complex region and the kayak- wagging region, with paramters 7 = 4.0 and A^ = 1.2. In 
part (a), the coupling constant e= 0.5, in part (b) e = 0.15 and in part (c) e = 0.1. All of 
these show qualitatively similar temporally intermittent behaviour. The fourier transform 
of the time series of stress calculated at a generic site and plotted on a doubly logarithmic 
scale is shown in Fig. (JTj). The spectrum in (c) of Fig. ((Tj) fits the relation P{f) ~ 1//^- 
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FIG. 4: [Color Online] Local dynamics in a ring of 200 lattice points, with coupling e = 0.1, showing 
the temporal evolution of si, S2 and n^. These are displayed for for (a) Xk = 0.9 and 7 = 2.0 (b) 
Afc = 1.1 and 7 = 2.0, and (c) X/. = 1.4 and 7 = 2.0. These points are drawn from the KT and KW 
part of the local phase diagram and represent states in which the director exhibits out-of-plane 
fluctuations i.e. Uz 7^ 0. However, si and S2 continue to exhibit regular, periodic oscillations. 

B. Spatio-temporal coherence and dynamics 

In order to quantify the degree of spatial coherence, we calculate the following quantity 
for the one- dimensional lattice: 
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We have calculate such a spatial coherence parameter for one specific component of the 
vector (ao, . . . , 04); however, qualitatively similar results are obtained for other components 
as well as for the full local stress, in the C region. When d tends to zero the degree of 
synchronization of the local variables is very high. On the other hand large d indicates low 
spatial synchronization, arising from a wide distribution of values of the local variables in the 
lattice. This quantity thus serves as a global order parameter characterizing the smoothness 
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FIG. 5: [Color Online] Time Period (T) vs. 7, for the case of a ring of 200 lattice points, with 
coupling e = 0.1 and with 7 varying across the T and KT regions of the phase diagram. The 
numerical results are represented as points and the fitted curve is: T ~ ^. 



of the spatial patterns exhibited by the evolution of the map. 

Fig. dH]) shows the time average of the deviation < d > oi a^ from the average value oq. 
To compute this, we first calculate the instantaneous deviation d via Eqs. ( 1T^ and (fT3|) . and 
then find the long-time average of this quantity. The spatial profile of the regular region with 
low < d > is characterised either by spatiotemporal fixed behaviour with all sites aligned, or 
spatial uniformity and temporal periodicity. There are also cases in the regular region where 
the sites, though not completely synchronized in space, are nevertheless phase synchronized. 

1. Spatio-temporal dynamics in the regular region 

Fig. ([9]) displays the space-time plot for 7 = 4, A^ = 1.1 and coupling constant e = 0.1 
(a) and 0.5 (b). The x-axis displays the lattice index and time is shown on y-axis, increasing 
from top to bottom. The profile is not spatially uniform and periodic in time for very weak 
coupling. As the coupling is increased, the system acquires spatial coherence and temporal 
periodicity. 
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FIG. 6: [Color Online] Local dynamics of a ring of 200 lattice points, with A^ = 1.2 and 7 = 4.0, 
showing the time evolution of si, S2 and Uz- These are shown for (a) e = 0.5 (b) e = 0.15 and 
(c) e = 0.1, illustrating behaviour in the complex or C regime. Note that regular time-periodic 
behaviour is favoured at large values of the spatial coupling constant e, following an initial transient. 

2. Spatio-temporal dynamics in the com,plex region 



The spatiotemporal behaviour of a representative case in the C or complex region is 
displayed in Fig. flTOj) . where 7 = 4, A^ = 1.17 and coupling constant e = 0.1 (a) and 0.5 
(b). It is evident that the space-time profile splits into bands, i.e. clusters of synchronized 
sites, where the local dynamics is either fixed (blue) or time-periodic (stripes). As we 
increase A^ (with 7 = 4) in Figs. ( HOl - [T2l) the length scale of the spatio-temporally 
intermittent pattern increases, finally yielding to the aligned region. This progression from 
frozen localized kinks/domains of fixed points in the spatial background of time-periodic 
behaviour, to infective bursts bearing the signature of spatiotemporal intermittency, is seen 
in many systems 53| , and often arises from a competition of fixed point patterns and time- 
periodic and quasi-periodic patterns. 
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FIG. 7: [Color Online] Log-Log plot of the Fourier transform vs. frequency, for (a) A^ = 1.17, (b) 
Afc = 1.20 and (c) Afc = 1.24. Here e = 0.1 and 7 = 4.0. The lattice is a ring of 200 sites. Note 
that for Afc = 1.17 (the wagging region), the regular oscillations show up as a delta function in the 
fourier transform. In the C or complex region, a smooth distribution of frequencies is seen, with a 
1/f falloff. 




FIG. 8: [Color Online] Average deviation from the mean value for e = 0.1 (a) and e = 0.5 (b), 
with Afc on the x-axis and 7 on the y axis. Note that large fluctuations (roughness) are seen in 
the KW and C regions. These data are for the 1-d system wiith the number of sites N = 200 and 
parameters as indicated on the figure. 
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FIG. 9: [Color Online] Density plot of order parameter part of shear stress. Here A^ = 1.1, 7 
=4.0 and e = 0.1 (a) and 0.5 (b). Time steps (running from top to bottom) are on the y-axis, 
and the lattice site index (i = 1,200) is on the x-axis. These figures represent space-uniform and 
time-periodic states, obtained using parameter values corresponding to the T region of the phase 
diagram. 















^(a)l 





\ - 


^ m 












?nnn 


- — 


■ ; 


■ 


- 




■ —- 


- 




'■ - ■ 










































4000 


: 


: ^E 




-■ 


6000 





■■= : 






8000 




. — - 








" ■ 


- — T~ 




























— m 









I 



50 



100 

L 



150 



I 



200 



2000 



4000- 



6000 



8000 



10000 




FIG. 10: [Color Online] Density plot of order parameter part of shear stress. Here A^ = 1.17, 7 
=4.0 and e = 0.1 (a) and 0.5 (b). Time steps (running from top to bottom) are on the y-axis, 
and the lattice site index (i = 1,200) is on the x-axis. These figures represent space non-uniform 
and time-periodic states, obtained using parameter values corresponding to the KW region of the 
phase diagram 
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FIG. 11: [Color Online] Density plot of order parameter part of shear stress. Here Afc = 1.20, 7 
=4.0 and e = 0.1 (left) and 0.5 (right). Time steps (running from top to bottom) are on the y-axis, 
and the lattice site index {i = 1,200) is on the x-axis. These figures illustrate how time-periodic 
regions are interspersed with domains of fixed point behaviour, reminiscent of spatiotemporal 
intermittency. The parameter values corresponding to the C region of the phase diagram, in a 
regime where the chaos is weak. 
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FIG. 12: [Color Online] Density plot of order parameter part of shear stress. Here Afc = 1.24, 7 
=4.0 and e = 0.1 (a) and 0.5 (b). Time steps (running from top to bottom) are on the y-axis, and 
the lattice site index (i = 1, 200) is on the x-axis. These figures illustrate non- uniform, time- varying 
states characteristic of spatio-temporally chaotic behaviour. The parameter values correspond to 
the C region of the phase diagram, in a regime where the chaos is strong. Note that larger values 
of e lead to larger and more uniform spatial structures. 
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FIG. 13: [Color Online] Density plot of order parameter part of shear stress. Here A^ = 1.24, 7 
=3.0 and e = 0.1 (a) and 0.5 (b). Time steps (running from top to bottom) are on the y-axis, and 
the lattice site index (i = 1, 200) is on the x-axis. 

V. THE TWO-DIMENSIONAL COUPLED MAP LATTICE 
A. Local temporal behaviour 

We have also explored the dynamics on a two-dimensional lattice. In the regular regions 
of the phase diagram, corresponding to the T,W, KT and KW states, the temporal behavior 
is very similar to that of the one dimensional case and is thus not shown separately. We 
thus concentrate on behavior in the complex or C region. 

Representative data showing the local temporal dynamics of the complex region is given 
in Fig ( 1T^ . They show chaotic behaviour, and there appears to be no qualitative difference 
between the one dimensional and two dimensional lattice cases. As in the one-dimensional 
lattice, increased coupling strengths suppress the chaotic region. The log-log plot of the 
fourier transform is shown in Fig. ( ITSll : a similar fit to P{f) ~ /^^ of the frequency spectrum 
of the stress can be obtained, as in the one-dimensional case. 
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FIG. 14: [Color Online] Temporal evolution of si, S2 and n^ in a two dimensional lattice of size 
102 X 102, with e = 0.1, for (a) Afc = 1.17 and 7 = 4.0 (b) A^ = 1.20 and 7 = 4.0 (c) A^ = 1.24 
and 7 = 4.0. All these state points are drawn from the C region of the local phase diagram. Note 
the existence of temporally intermittent behavior, analysed in terms of its frequency spectrum in 

Fig.m 

B. Spatio-temporal behaviour 

To quantify the degree of spatial coherence in 2- dimensional lattices we calculate the 
quantity, generalizing from the one-dimensional case studied in an earlier section: 
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Again, as in the 1-dimensional case, when d tends to zero the degree of synchronization of the 
local variables is very high. On the other hand large d indicates low spatial synchronization, 
and arises from a wide distribution of values of the local variables in the lattice. 
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FIG. 15: [Color Online] Log-Log plot of the absolute value of the Fourier transform of the stress 
vs. frequency uj, for (a) Afc = 1.17, (b) A^ = 1.20 and (c) Xk = 1.24. Here e = 0.1 and 7 = 4.0. The 
lattice contains 100 x 100 sites. Note the relatively smooth background, indicating the presence of 
a continuous set of frequencies. The fall-off is consistent with a l/o;^ behaviour. 

Fig. flT6l) shows the space-time average of the deviation defined in Eqn. flTSj) . The left 
panel (a) displays results for e = 0.1 and the right panel (b) for e = 0.5. It is clear that 
higher coupling strengths make the system more uniform in space. Also, it appears that the 
regions with kayak-tumbling, kayak-wagging and complex local dynamical behaviour show 
more deviation in the spatial profile, exhibiting more spatial inhomogeneity. 

1. Regular regime 

In this section we discuss the spatial profile of our coupled map lattice in two dimensions. 
As in the one- dimensional case, we start with random initial conditions and analyze the 
space profile after omitting a transient regime. We analyse the density plots of the shear 
stress contribution to the order parameter, in different dynamical regions. Considering the 
space-time behaviour of the system in the regular region, with local dynamics belonging 
to the aligned, wagging, tumbling and kayak-tumbling region, reveals spatially uniformity 
states which are periodicity in time. These are closely related to the states obtained in the 
oine-dimensional case and are not discussed further here, as we will concentrate on results 
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FIG. 16: [Color Online] Average deviation d (see text) from the mean value of a^ for e = O.l(riglit) 
and e = 0.5(left). The quantity A^ is plotted on the x-axis and 7on the y-axis. The lattice is a 
100 X 100-site lattice. Note that this roughness is largest in the KT, KW and C regions, especially 
for large values of 7. As e is increased, the roughness decreases, as increasing spatial homogeneity 
is promoted. 

obtained in the physically more interesting C regime. 



2. Complex regime 

The configurations in Figs. (fT7|) and (fT8|) are from the complex region. When the coupling 
becomes very large, one obtains spatially uniform states, as clear from Figs. flT7|) and( [T8l) . 
In Fig. [191 we have chosen points (a) and (b) from the KT region of the local phase diagram 
and (c) from the C dynamical region. After leaving 10'^ transient steps , we have plotted 
one row of a 100 x 100 lattice at a single time instant. On the x axis we plot e and on 
the y axis we plot the stress at 100 points of the lattice at one time step. It is evident 
that for high coupling strength e, the system goes to a space-synchronized state. For low 
coupling constants, on the other hand, there is a typically wide distribution of stress values 
at different sites, indicating spatial inhomogeneity. 
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FIG. 17: [Color Online] (a) A color plot of order parameter part of shear stress. The quantities A^ 
= 1.25, 7 =4.0 and e = 0.1. The time t is plotted on the y-axis, which depicts the time evolution 
of the stress computed on one row (x-axis) of the 100 x 100 lattice, (b) shows a snap shot of the 
full lattice at an intermediate time step. 

VI. QUANTIFYING SPATIO-TEMPORAL COMPLEXITY 



In this section, we report results which quantify the spatio-temporal complexity in the 
one-dimensional coupled map lattice specified in Eq. [61 To understand the nature of the 
complex behaviour represented in the phase diagram, we perform calculations of the spec- 
trum of Lyapunov exponents, as shown in Fig. [201 We first choose several values of the 
parameters A^ and 7 within the complex region and evolve the coupled map. After waiting 
for an initial number of time steps to eliminate transients, we calculate the Jacobian matrix 
at each time step. We then consider a small deviation from the attractor and iteratively 
multiply this deviation by the Jacobian, orthonormalizing this vector at each time step. 
From this we calculate the Lyapunov exponent, using the method described in Ref. |54l ]. 

These results are illustrated in Fig. [201 which exhibits the values of the first four Lyapunov 
exponents, computed for parameter values A^ and 7 for system sizes L = 3, 10, 50 and 100, 
as a function of the coupling constant e. As can be seen from these figures, our results are 
the following. Qualitatively, in the complex regime, the first Lyapunov exponent is always 
positive, even as the system size and the spatial coupling are increased. The local value 
of this exponent is also positive. This value decreases further with spatial coupling but 
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FIG. 18: [Color Online] (a) A color plot of order parameter part of shear stress. The quantities A^ 
= 1.25, 7 =4.0 and e = 0.5. The time t is plotted on the y-axis, which depicts the time evolution 
of the stress computed on one row (x-axis) of the 100 x 100 lattice, (b) shows a snap shot of the 
full lattice at an intermediate time step. Note that increased spatial coupling favours the space 
uniform time-periodic state. 
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FIG. 19: [Color Online] Value of stress of one row of the lattice at one instant of time on the y 
axis as a function of the coupling strength e on the x axis. Note that the broad spectrum of local 
stress values seen at small values of the coupling contract to an essentially unique value at large e. 
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FIG. 20: [Color Online] Plots of the Lyapunov exponents obtained from our calculation (see text) 
for different system sizes as a function of the coupling constant. The figures represent (a) he first 
Lyapunov exponent Ai (b) the second Lyapunov exponent A2 c) the third Lyapunov exponent A3 
d) the fourth Lyapunov exponent A4. Note that all these exponents tend to zero in the limit of 
large system sizes and that increasing the coupling between sites has the effect of reducing the 
Lyapunov exponent if it were initially positive at zero coupling and increasing it, if it were initially 
negative. 
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remains positive. Roughly speaking, larger lattice sizes show larger values for this exponent, 
consistent with results from the one-dimensional PDE calculation. These values appear to 
saturate for small coupling values but decrease for larger values of the spatial coupling. 

The second and higher order Lyapunov exponents, in our calculation, are small and 
negative for the smallest lattice sizes, but move to values that are close to zero as the lattice 
size is increased. At small couplings, for the larger lattices, this value is positive but goes 
negative as the coupling strength is increased. Thus, the data for the Lyapunov exponents 
are consistent with the general conclusion that going to larger lattice sizes stabilizes chaos, 
whereas increasing the coupling between sites suppresses complex spatial behaviour. The 
clustering of Lyapunov exponents around zero in the large system size limit is consistent 



with the emergence of spatio-temporal intermittency on large scales 55N57|. 



VII. CONCLUSIONS 

In summary, this paper reports a study of a coupled map lattice model constructed to 
study spatio-temporal aspects of rheological chaos in sheared nematic solutions. Our study 
was based on the construction of a suitable local map capable of reproducing the physics of 
the spatially uncoupled (equivalently, uniform) limit, including the variety of phases and the 
complex phase diagram obtained for that case. When such maps are placed on a regular (one 
or two-dimensional) lattice and coupled diffusively through a variety of coupling schemes, 
the local dynamics in the coupled map shares close similarities with the map describing the 
spatially uncoupled case. 

Our studies of the coupled map in both one and two dimensions indicates that regimes of 
regular behaviour largely exhibit space-uniform and time-periodic states, with the coupled 
dynamics roughly following the uncoupled case. We have analysed the dynamical behaviour 
of the two quantities which characterize local order in the nematic, the uniaxial Si and the 
biaxial S2 order parameters, examining their time evolution in the different states. 

In contrast, in the complex or C region of the local phase diagram, such coupling leads 
to states that exhibit spatio-temporal intermittency and chaos. We have characterized such 
states by examining the Lyapunov spectra as well as the frequency dependence of the time 
series of physical quantities such as the stress. We find evidence for a broad, power-law 
distribution of time-scales in the problem. Further, in the complex region, one often sees 
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a coexistence of regular (lamina) and chaotic regimes as a prelude to fully developed chaos 
in which dynamical fluctuations occur independently from site to site. In some regimes, 
periodic bands immersed in a more complex, fluctuating background are obtained, suggestive 
of the possibility of transient shear bands stabilized by the dynamics, a feature also present 



in ODE-based studies of this problem 27|, |28|, |35| . The basic scale of these complex dynamical 



patterns is alterable by changing the coupling constant, indicative of self- similarity in the 
spatio-temporally intermittent case. At very large values of the coupling constant, the space 
profile is expected to become uniform; however, for small and intermediate values of this 
coupling constant, the spectrum of Lyapunov exponents merges to zero, consistent with our 
observation of generic spatio-temporal intermittency in the weak coupling case. 

We have experimented with using spatial coupling terms which represent the advective 
effects of the shear flow, coupled to flxed boundary conditions where the orientation and 
magnitude of the order parameter are flxed at the boundary. Such terms appear, at small 
amplitude, to mainly distort the sorts of dynamical structures obtained for the symmetric 
coupling state and seem to evolve smoothly from them. 

The usefulness of coupled map lattice representations of the spatio-temporal dynamics of 
systems exhibiting chaos in their local dynamics is that such representations often provide 
both useful physical insights as well as are computationally easier to simulate than their 
PDE versions. In that sense, the problem of rheochaos in sheared nematics offers an ideal 
setting for CML methods, since the local dynamics of the sheared nematic is highly non- 
trivial, exhibiting a variety of temporally periodic as well as chaotic states. As shown here, 
the variety of non-trivial spatio-temporal behaviour exhibited by sheared nematics is very 
largely a consequence of simply coupling these dynamical degrees of freedom in space. The 
physics appears substantially independent of how precisely this spatial coupling is done, 
with the simple lattice model with parallel update exhibiting virtually all the behaviour of 
the more complex and computationally intensive studies of the appropriate PDE's. This, 
together with the speciflc results presented in this paper for our coupled map approach to 
rheochaos in sheared nematics, is our central conclusion. 

Further, order-parameter-based models, such as the one described in this paper and in 



the work of Refs. 



27 



28 



35| . contain essential non-linear terms in the free energy. It is 



these terms that are responsible for the non-trivial local dynamics captured in our local 



map as well as in the coupled map lattice. Ref. 35| emphasizes the role of "additional 

29 



complex collective dynamics" arising from such nonlinearities which is not captured in the 
DLS model but is relevant to the qualitative nature of the intermittent and chaotic behaviour 
seen in this system. Such non-linearities are naturally accounted for in our approach. 

Our study of the spatio-temporal dynamics of sheared nematics using CML methods pos- 
sibly represents the first extension of such methods to the problem of rheochaos. In contrast 
to previous work based on ODE's which studied only the one-dimensional case, it is relatively 
easy to extend our CML methodology to higher dimensions, even to the experimentally rel- 
evant three-dimensional case. It would be interesting to see how, if at all, hydrodynamic 
effects can be incorporated in models of this type. Whether other experimental systems of 
sheared complex fluids which exhibiting rheochaos can be fruitfully analysed using similar 
coupled map approaches remains to be seen. 
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